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The family of the superconducting quasi-skutterudites (Caa;Sri_ 2 ,) 3 Rh 4 Sni 3 features a structural 
quantum critical point at Xc = 0.9, around which a dome-shaped variation of the superconducting 
transition temperature Tc is found. Using specific heat, we probe the normal and the supercon¬ 
ducting states of the entire series straddling the quantum critical point. Our analysis indicates a 
significant lowering of the effective Debye temperature on approaching Xc, which we interpret as a 
result of phonon softening accompanying the structural instability. Furthermore, a remarkably large 
enhancement of 2 A/A:_bTc and AC/yTc beyond the Bardeen-Cooper-Schrleffer (BCS) values is found 
in the vicinity of the structural quantum critical point. The phase diagram of (Caa;Sri_a;) 3 Rh 4 Sni 3 
thus provides a model system to study the interplay between structural quantum criticality and 
strong electron-phonon coupling superconductivity. 

PACS numbers: 65.40.Ba, 74.40.Kb, 74.25.-q, 74.62.-c 


The interplay between superconductivity and struc¬ 
tural instability has been an important theme in con¬ 
densed matter physics. The proximity to structural in¬ 
stability has been suggested to play a role in several su¬ 
perconductors with relatively high transition tempera¬ 
ture (Tc) HHZIj with recent examples including IrTe 2 i- 
[TO] as well as Fe- and Ni-based superconductors [TTHIK] . 
To examine the prospect of stabilising or even enhancing 
superconductivity near structural instability, access to 
superconducting systems featuring a second-order struc¬ 
tural phase transition, which can be tuned away, is highly 
desirable. 

A large family of superconducting stannides [ig [i2] 
with composition AsTqSnia, where A=La,Sr,Ca and 
T=Rh,Ir, has recently attracted considerable attention 
[T5H55] . Tc in these compounds are typically less than 
10 K, and the superconducting gap function is established 
to be nodeless at the Fermi surface [I9H22] . What makes 
these systems interesting is the discovery of an additional 
second-order structural phase transition at T* in some of 
the compositions [5DH5^ . Systematic variation of T* has 
been observed when the unit cell volume of the crystal 
is varied via chemical substitution or applied pressure, 
resulting in phase diagrams suggestive of a strong inter¬ 
play between the structural and superconducting orders, 
which is in striking resemblance to the phase diagrams 
of heavily studied superconductors found on the border 
of magnetism . 
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SralrqSnis and Sr 3 Rh 4 Sni 3 are two prominent mem¬ 
bers of the family which exhibit both T* and Tc- 
Clear lambda-like jump in the specific heat can be seen 
at T* = 147 K and 138 K in Sr 3 lr 4 Sni 3 |3D] and 
Sr3Rh4Sni3|32j. respectively. Crucially, T* can be sup¬ 
pressed rapidly by applying pressure, or substituting cal¬ 
cium by strontium [an Eg. In the Ir substitution series, 
i.e. (Caa:Sri_ 3 ;) 3 lr 4 Sni 3 , T* reaches as low as 33 K when 
X = 1. Full suppression of T* can only be achieved by 
applying 18 kbar to Ca3lr4Sni3[3T]. In the Rh substitu¬ 
tion series, T* can be fully suppressed solely with calcium 
substitution, giving rise to a structural quantum critical 
point (QCP) at ambient pressure at a: « 0.9 [32]. Con¬ 
current with the suppression of T* , Tc rises gently in both 
Ir- and Rh-series, and reaches the maximum value near 
where T* —> 0 Ell El] ■ 

A second-order structural phase transition involves the 
softening of the phonon mode that is responsible for the 
structural order [40l |4T| . This will inevitably affect the 
low energy excitations of the system. In addition, the 
nodeless nature of the superconducting gap function sug¬ 
gests a conventional Cooper pairing mechanism medi¬ 
ated by phonons, allowing the discussion of the coupling 
strength. The tunability of the (Ca 4 ;Sri_a;) 3 Rh 4 Sni 3 se¬ 
ries across the structural QCP without applying high 
pressure offers an important opportunity to understand 
the underlying lattice dynamics and the relationship be¬ 
tween structural and superconducting orders. In this Let¬ 
ter, we report the results of our heat capacity measure¬ 
ments on (Caa;Sri_a;) 3 Rh 4 Sni 3 for seven Ca concentra¬ 
tions, spanning the range from x = 0 to a: = 1. The 
data enable us to extract important parameters rele- 
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FIG. 1. (Color online) The structure of Sr 3 Rli 4 Sni 3 (a) above 
and (b) below T*. The Sr, Rh, and Sn atoms are represented 
by spheres with descending sizes, and icosahedra, nntilted and 
tilted trigonal prisms are colonred in tangerine, blue and pink. 
Crystallographically distinct Sn(2) atoms are coloured green, 
red, yellow and orange, and a single unit cell of each structure 
is indicated, (c) and (d) show the subtle tilting of the trigonal 
prisms RhSn(2)6 around an Sn(l)Sn( 2 )i 2 icosahedron. 

vant to both superconducting and structural phase tran¬ 
sitions, which reveal a remarkable trend on approaching 
the structural QCP. 

Single crystals of (Caa:Sri_ 2 ,) 3 Rh 4 Sni 3 were synthe¬ 
sized by Sn flux method using similar parameters as de¬ 
scribed in ref. |18| . Heat capacity was measured in a 
Physical Property Measurement System (Quantum De¬ 
sign) using a standard pulse relaxation method. The 
background heat capacity of the sample platform and 
Apiezon N-grease, which was used as adhesive, was care¬ 
fully measured from room temperature down to the low¬ 
est attainable temperature. To accurately determine the 
normal state heat capacity of the samples, additional 
background measurements were performed at high mag¬ 
netic field (< 5 T). The mass of the samples ranges from 
^5 mg to ^50 mg. 

Above T*, Sr 3 Rh 4 Sni 3 adopts a cubic Pm3n structure 
(Fig. 1(a)) featuring Sn(l)Sn( 2 )i 2 icosahedra connected 
via RhSn(2)6 trigonal prisms, with the Sr cations occu¬ 
pying the space between icosahedra. In this high tem¬ 
perature centrosymmetric phase, the bond lengths in the 
Sn(l)Sn( 2 )i 2 icosahedra are identical. Below T*, these 
icosahedra distort leading to each having four groups of 
Sn(l)-Sn(2) bonds of different lengths |32]. This requires 
six of the eight trigonal prisms around an icosahedra to 
tilt, with only those along the [ 111 ] direction remain¬ 
ing untilted. This results in a symmetry lowering to a 
non-centrosymmetric IA3d cubic structure with a dou- 



FIG. 2. (Color online) (Left) The specific heat of Sr 3 Rh 4 Sni 3 
{x = 0 ) and Ca 3 Rh 4 Sni 3 (a: = 1 ) as a function of temperature 
at zero external magnetic field (closed symbols). T* indicates 
the temperature at which a structural transition takes place. 
(Right) A closeup of the specific heat at low temperatures. 
The midpoint of the jump defines Tc, which can be suppressed 
by applying external magnetic field (open symbols). 

bled unit cell (see Fig. 1(b)). 

Fig. i shows the temperature dependence of the 
specific heat for Sr 3 Rh 4 Sni 3 and Ca 3 Rh 4 Sni 3 . In 
Sr 3 Rh 4 Sni 3 , a lambda-like jump in the specific heat oc¬ 
curs at T*~138 K, which corresponds to a second-order 
structural phase transition [32]. At low temperatures, 
another jump in the specific heat is detected (see the 
right panel of Fig. [^, which marks Tc and is in excellent 
agreement with previous electrical resistivity measure¬ 
ment |32| (see Fig. SI of Supplemental Material). In 
Ca 3 Rh 4 Sni 3 , the specific heat varies smoothly except at 
Tc = 8 K. This is again consistent with electrical resis¬ 
tivity measurement which did not detect any evidence of 
T* (Fig. SI). 

The specific heat jump at Tc can be suppressed by 
magnetic field, thereby exposing the normal state. The 
low temperature specific heat in the normal state can be 
adequately described by the following model: 

„0e/T / 0 p \ ^ 

C„(T) = 7T + /3T3 + «^-^-^^-^(^^j , (1) 

where the first term is the contribution from conduc¬ 
tion electrons, the second and the third terms are due 
to phonons with dispersion relations uj (x q (Debye) and 
q-independent uj (Einstein), respectively. Qe is the Ein¬ 
stein temperature jUj. In the inset to Fig. |^a) where 
Cn/T is plotted against T^, the lowest temperature data 
follow a straight line and the value of 7 and /3 could be ex¬ 
tracted. The excess specific heat at higher temperatures 
is attributed to the Einstein term. Using this model, we 
successfully described the low temperature specific heat 
of all compositions studied (main panel of Eig. j^a) and 
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FIG. 3. (Color online) (a) Temperature dependence of the 
specific heat of (CaiSri_a;) 3 Rh 4 Sni 3 in the normal state for 
X = 0,0.25,0.5,0.75,1. Inset shows C/T as a function of . 
The experimental data (open symbols) at low temperatnres 
are fitted linearly (solid lines). 7 and j3 in Eq. 0 are ex¬ 
tracted from the intercept and slope of the fit, respectively. In 
the main panel, the measnred specific heat is fitted nsing the 
full expression of Eq. 0 with the extracted values of 7 and f5. 
(b) The specific heat of the critical composition {x = 0.9) in 
zero (closed symbols) and an applied (open symbols) external 
magnetic field as a function of temperature. The solid curve 
shows the fitting of the data using Eq. 0 and the value of 
the superconducting gap can be estimated. Inset shows the 
plot of C/T as a function of temperature near T^. The data 
are linearly extrapolated around Tc to determine the jump in 
the specific heat AC/Tc = {Cs — Cn)/Tc. (c) The rescaled 
electronic specific heat in the superconducting state, plotted 
on the logarithmic scale, as a function of T~^ for x = 0 and 
0.95. The straight lines from T^/T = 0.8“^ = 1.25 are guides 
to the eyes, with slopes a related to A via a = A/ksTc- 


Fig. S6-S7 in Supplemental Material). 

Next, we extract the parameters relevant to supercon¬ 
ductivity. We take the midpoint of the specific heat 
anomalies at low temperatures as T^. Following the 
scheme shown in the inset to Fig. [^b), we estimated the 
normalized specific heat jump AC/^T^ with the value 
of 7 obtained from the analysis of Cn described above. 
In addition, we estimate the superconducting gap A by 
analysing the zero field heat capacity data well below 


using 


Cs (T) = A exp 



A 


+/3r3-bAc. 


gOE/T 

(gOETr A T)2 



( 2 ) 

The phonon contributions from the the second and the 
third terms in Cg are assumed to be identical to that in 
Cn- Hence, A and A are the only adjustable parameters 
in the fitting. The usage of a single s-wave gap is justified 
based on previous studies, in which no evidence of low 
energy quasiparticle excitation was detected in the super¬ 
conducting state, indicating a nodeless superconducting 
gap function [IMH EH- For T ^ Tc, the BCS the¬ 
ory predicts that the gap closes in a fashion ~ \/Tc — T 
[43l m]. Therefore, we fit the zero field heat capacity 
data below O.STc with Eq. <0 for all seven samples (see 
Fig. [^b) for X = 0.9, and Fig. S2-S8 for others). The 
validity of our model is clearly demonstrated in Fig. |^c) 
where the electronic part of the specific heat at zero field, 
Ces, is plotted on the logarithmic scale as a function of 
T“^, with the solid lines drawn using the values of A 
extracted from the analysis of CHTipS]. 

The results of these analyses are summarised in Fig. 
4. In Fig. 4(a), we overlay T* and Tc determined from 
the current heat capacity measurement on the universal 
phase diagram constructed earlier [32]. For x > 0.5, the 
T* anomaly becomes significantly smeared out in specific 
heat, although it is still visible in electrical resistivity. For 
all compositions studied here, Tc values are in excellent 
agreement with previous reports [H] [32] . The phase di¬ 
agram constructed clearly shows the suppression of T* 
with increasing Ca content, giving rise to a structural 
QCP at a: = 0.9. A broad dome-like variation in Tc is 
found to peak around the QCP. 

We adopt the viewpoint that the second-order struc¬ 
tural phase transition is accompanied by the softening of 
the relevant phonon mode frequency uJs,q on approach¬ 
ing T* from T > T*, and hence the branch of uj{q) on 
which this mode lies will be affected. In Landau’s theory, 
Wg Q oc {T — T*) plays the role of the inverse generalised 
susceptibility [40] |4T] . Below T *, the structure relaxes to 
a new stable configuration, and (jJs^q will harden again. 
This temperature dependence of the soft phonon mode 
frequency across a structural transition has been most 
clearly observed in ferroelectric compounds, e.g. SrTiOa 
pnHKT] . For the present system, T* = 0 K at the QCP. 
Consequently, at the QCP, uJs,q —)■ 0 at 0 K and we are 
left with an abundance of low-lying phonon modes at low 
temperatures. 

First evidence of the existence of the soft mode comes 
from the electrical resistivity measurements of both 
(Caa;Sri_ 3 ;) 3 Rh 4 Sni 3 and (Caa;Sri_ 3 ;) 3 lr 4 Sni 3 . At the 
structural QCP, namely in (Cao, 9 Sro.i) 3 R-h 4 Sni 3 at am¬ 
bient pressure and Ca 3 lr 4 Sni 3 under 18 kbar, a distinct 
linear-in-T variation of the electrical resistivity was ob¬ 
served over a wide temperature range, which was at- 
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FIG. 4. (Color online) (a) The phase diagram of 

(Caa;Sri_a;) 3 Rh 4 Sni 3 on the T — x plane. Polygon symbols 
represent data extracted from heat capacity measurement. 
Cross symbols are data from electrical resistivity measure¬ 
ment in Ref. [32]. The region beyond a; = 1 is accessible 
by applying pressure. Right panels show the x dependence 
of (b) the effective Debye temperature and the Einstein tem¬ 
perature, (c) 7 , (d) ACj'^Tc and 2A/fesTc, respectively. The 
numerics of these parameters are tabulated in Table SI of the 
Supplemental Material. 


tributed to a strong coupling to the soft phonon modes 
[3III32]. Lattice dynamics calculations further revealed 
that the soft mode is located at M(0.5, 0.5, 0) 15^ . 

which is consistent with the structural refinement per¬ 
formed using data from x-ray diffraction studies |3T1 [32] , 
and a recent inelastic scattering results in Ca^IrASnmfS^]. 
In Ca 3 Rh 4 Sni 3 , which is beyond the structural QCP and 
hence does not undergo a structural transition, the soft 
mode frequency is 0Js,q 0.35 THz [32]. Since uJs,q is a 
local minimum at M, we can Taylor-expand around 

and obtain to oc ^+ to^ where q' is the wavevector 
measured from M. We emphasise that the lattice dynam¬ 
ics calculations were performed for 0 K. Thus, to induce a 
structural quantum phase transition from the quantum 
disorder side, one needs to reduce uJs,q to zero, result¬ 
ing in a phonon dispersion relation to ot q'. Hence, at the 
QCP, an additional, acoustic-like phonon branch emerges 
at M. 

We are now in the position to explain the x-dependence 
of the remaining quantities as shown in Fig. and to put 
the results in the context of structural quantum critical¬ 
ity associated with this system. First of all, the contri¬ 
bution of the additional acoustic-like branch to the spe¬ 
cific heat will be proportional to following Debye’s 
treatment. Indeed, we found that /3 (c.f. Eq. (1)) is 


significantly larger, or equivalently, the effective Debye 
temperature is considerably suppressed, at around 
X = 0.9 [54]. Our analysis indicates a rapid decrease in 
0|f as X increases, with a broad minimum around x = 0.9 
where 0|f is ~ 33% lower than that at x = 0. Con¬ 
current with the strong variation of 0|^, a qualitatively 
similar variation is observed in Oe- This is consistent 
with our expectation that an abundance of low-energy 
phonon modes become readily populated as a result of 
the softening of tOg^Q. 

From Fig. 4(d), we see that in addition to an enhance¬ 
ment in Tc with increasing x, both /S.C and 2A/kBTc 
are appreciably larger on the Ca-rich end of the phase 
diagram. Criticality in ACjjTi. can be clearly observed 
near the structural QCP at x = 0.95. Except for x = 0, 
ACj'^Tc is noticeably larger than the BCS value of 1.43 
[431 l44l 133] for all samples studied, and reaches a value 
as high as ^3.47 at x = 0.95. Note that both Tc and 
7 increase with increasing x (Figs. |^a) and (c)), there¬ 
fore the enhancement in AC/qTc is accompanied by a 
much larger AC. Moreover, 2A/kBTc also experiences 
a sizeable enhancement near the structural QCP, reach¬ 
ing a remarkably high value of ~6.67 at x = 0.95. In 
fact, theory of strong coupling superconductivity give 
|AC/ 7 T,U„, ~ 2.98 and \2A/kBTc\max ~ 5.57 [33], 
which are both smaller than our observed values near the 
structural QCP. Similar observation has recently been re¬ 
ported from /rSR studies, which found a large 2A/kBTc 
of ~ 8 near the pressure-induced structural quantum crit¬ 
ical point in the sister compound Ca 3 lr 4 Sni 3 [33]. 

Combining all parameters relevant to the normal and 
superconducting states from systematic studies across 
the phase diagram of (Caa;Sri_ 2 ;) 3 Rh 4 Sni 3 , an intrigu¬ 
ing picture emerges: at x = 0, we have a superconductor 
close to the weak coupling limit, with both 2A/kBTc and 
ACfjTc close to the expected ratios in the BCS theory 
[431|44l[33]. This is in stark contrast to the Ca-rich end, 
in the vicinity of the structural QCP, where these ratios 
are signihcantly enhanced. Thus, the phase diagram can 
be regarded as the outcome of fine tuning the coupling 
strength, which ultimately give rise to the intricate in¬ 
terplay between the phases observed. 

In summary, we have analysed the heat capacity data 
of (Caa:Sri_ 2 ,) 3 Rh 4 Sni 3 , in the normal and supercon¬ 
ducting states, for seven x values straddling across the 
structural QCP at x = 0.9. The effective Debye tem¬ 
perature and the Einstein temperature are found to be 
significantly reduced near x = 0.9, which we attributed 
to the softening of a phonon mode responsible for the 
structural transition. Additionally, Tc is found to peak 
near x = 0.9, accompanied by a remarkable enhance¬ 
ment in 2A/kBTc and ACj^Tc far beyond the BCS val¬ 
ues. (Caa;Sri_a;) 3 Rh 4 Sni 3 thus serves as an unprece¬ 
dented model system for studying the intricate interplay 
between structural instability and strong coupling super¬ 
conductivity. 
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